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Abstract 

A phase field model is presented to study dislocation formation (coherency loss) in two-phase 
binary alloys. In our model the elastic energy density is a periodic function of the shear and 
tetragonal strains, which allows multiple formation of dislocations. The composition is coupled to 
the elastic field twofold via lattice misfit and via composition-dependence of the elastic moduli. By 
numerically integrating the dynamic equations in two dimensions, we find that dislocations appear 
in pairs in the interface region and grow into slips. One end of each slip glides preferentially into 
the softer region, while the other end remains trapped at the interface. Under uniaxial stretching 
at deep quenching, slips appear in the softer region and do not penetrate into the harder domains, 
giving rise to a gradual increase of the stress with increasing applied strain in plastic flow. 
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I. INTRODUCTION 



In crystalline solids many kinds of phase transformations are strongly influenced by the 
elastic Held^^^^. Since the first work by Cahn^ most theoretical studies have been focused 
on the coherent case in which the lattice planes are continuous through the interfaces. In 
the incoherent case, on the other hand, dislocations appear around the interfaces and the 
continuity is lost partially or even completely. Such incoherent microstructures emerge in 
various alloys when the lattice constants or the crystalline structures of the two phases are 
not close^iSiSiiS. Moreover, they are produced in plastic flow because dislocations generated 
tend to be trapped at the interfaces. In particular, coherency loss has been extensivesly 
studied in the presence of 7' precipitates of the LI2 structure2J£. 

Theory for the incoherent case is much more difficult than for the coherent case, obviously 
because the effects cannot be adequately described within the usual linear elasticity theoryii. 
The aim of this paper is hence to present a simple mathematical model reasonably describing 
the incoherent effects in binary alloys. Use will be made of a recent nonlinear elasticity theory 
of plastic flow by one of the present authorsi^. 

A number of authors have studied composition changes around dislocations fixed in space 
and time, which lead to a compositional Cottrell atmospheroi^ or preferential nucleation 
around a dislocation^^'—. As recent numerical work in two dimensions, phase separation has 
been studied by Leonard and DesaiiS^ and by Hu and Chenii using a continuum Ginzburg- 
Landau or phase field model in the presence of fixed dislocations. In these papers, dislo- 
cations preexist as singular objects before composition changes. We also mention atom- 
istic simulations of dislocation motion infiuenced by diffusing solutesi^ or by precipitated 
domainsi^. 

Mechanical properties of two-phase solids are very different from those of one-phase 
gQjj^g7ji3j2o_ ^Yie presence of precipitated domains, dislocations can be pinned at the 
interface regions and networks of high-density dislocations can be formed preferentially in 
softer regions after deformations^. These effects are very complex but important in technol- 
ogy. Our simulations will give some insights on the behavior of dislocations in two-phase 
states. 

This paper is organized as follows. In Section 2 we will present the free energy functional 
for the composition and the elastic field, in which the elastic energy density is a periodic 



function of the tetragonal and shear strains and the composition is coupled to the elastic 
field. In Section 3 we will construct dynamical equations. In Section 4 numerical results 
will be given on the dislocation formation around domains and on the stress-strain relations 
under uniaxial stretching. 



II. FREE ENERGY FUNCTIONAL 

We consider a binary alloy consisting two components, A and B, neglecting vacancies 
and interstitials. The compositions, ca and cb, of the two components satisfy ca + cb = 1- 
In real metallic alloys undergoing a phase transition, there can be a change in the atomic 
configuration within unit cells as well as in the overall composition, resulting in ordered 
domains with the so-called LIq or LI2 structure^'^'^^. However, in this paper, the composition 
difference is the sole order parameter, 

tp = ca- Cb, (2.1) 

for simplicity. The other variables representing the order-disorder phase transition are ne- 
glected. Then ijj is in the range —1 <1 and 

CA = + CB = lil-^). (2.2) 

In our free energy F = J drf the order parameter ip and the elastic displacement vector 
u = {ux, Uy) are coupled. The free energy density / is of the form, 

/ = /bw(V') + ^1^1' + "eiV- + /el. (2.3) 
The first term is the Bragg- Williams free energy density expressed as^ 

'°./BW = ^ln(l + ^) + i^ln(l-^) 



fceT^"" 2 

- ToV^V2T, (2.4) 

where Vq is the volume of a unit cell representing the atomic volume, Tq is the mean-field 
critical temperature in the absence of the coupling to the elastic field. If -C 1, we 
obtain the Landau expansion v^fsw /k-oT = (1 — To/T)'?/;^/2 + -?/;^/24 + ■ ■ ■ . However, we will 
not use this expansion form because we are interested in the deeply quenched case. The 
second term in (2.3) is the gradient term where C is a positive constant. The parameter 



a represents the strength of the couphng between the composition and the dilation strain 
ei = V ■ n. This couphng arises in the presence of a difference in the atomic sizes of the two 
species and is consistent with the empirical fact that the lattice constant changes linearly as 
a function of the average composition in many one-phase alloys (Vegard law). It gives rise 
to a difference in the lattice constants of the two phases in phase separation (lattice misfit). 
It also explains a composition inhomogeneity (Cottrell atmosphere in one-phase states or 
precipitate in two-phase states) around a dislocation. 
In two dimensions /ei depends on the following strains, 

ei = Vx.Wx + Vj^tiy, 

63 = VyUx Vx%, (2.5) 

where Vx = djdx and Vy = d/dy. The elastic displacement u is measured in a reference 
one-phase state at the critical composition. We call 62 the tetragonal strain and 63 the shear 
strain. In this paper we use a nonlinear elastic energy density of the form, 

/el = ^i^6? + $(^,62,63). (2.6) 

The first term represents the elastic energy due to dilation with K being the bulk modulus. 
The second term arises from anisotropic shear deformations defined for arbitrary values of 
62 and 63. Assuming a square lattice structurei^-, we set 

$ = 0[1 - cos(27r62)] + 0[1 - cos(27r63)]. (2.7) 

The principal crystal axes are along or make angles of ±7r/4 with respect to the x or y axis. 
In Fig.l we plot $ as a function of 62 and 63 for the case /i2 = /^s = /^o iii units of fiQ. If the 
system is homogeneous, elastic stability is attained for d'^^/de'^ > and d'^^/de^ > or in 
the regions \e2 — n\ < 1/4 and [63 — m| < 1/4 with n and m being integer values^^. 
For small strains |62| ^ 1 and I63I <^ 1, it follows the usual standard forn*ii, 

^ = lf^2el + ^fisel (2.8) 

in the linear elasticity theory. Therefore, 

/i2 = ^(Cll-Ci2), fX3 = Cu, (2.9) 



in terms of the usual elastic moduli Cu, C12, and C^^. In the original theorj^ the isotropic 
linear elasticity with constant /i2 = /^s was assumed. Subsequent theories treated the case 
of the cubic linear elasticity with constant /i2 and /x ^^i^i'^i^i^^ . In the present paper, while K 
is a constant, fi2 and fi^ depend on the composition as 

Ai2 = /^2o + Ai2i^, = /"so + Ai3i^- (2.10) 

If yU2i > and /isi > 0, the regions with larger (smaller) ip are harder (softer) than those with 
smaller (larger) ip. It is known that this elastic inhomogeneity gives rise to asymmetric elastic 
deformations in two-phase structures and eventual pinning of the domain growth^i^^^. 

In our theory 62, 63) in (2.7) is the simplest periodic function of 62 and 63 with period 
1. The periodicity arises from the fact that the square lattice is invariant with respect to a 
slip of the crystal structure by a unit lattice constant along a line parallel to the x or y axis. 
Notice that, under rotation of the reference frame by 6, 62 and 63 are changed to 65 and 63, 
respectively, within 

63 = 62 cos 26^ + 63 sin 2^, 

= -62 sin 2^ + 63 cos 2^. (2.11) 

For 9 = 7r/2 we have e'g = —62 and 63 = —63, so /ei in (2.7) remains invariant. For 6 = 7r/4 
we have e'g = 63 and 63 = —62 and recognize that the roles of tetragonal and shear strains are 
exchanged. For /i2 = /i3, the linear elasticity in (2.8) becomes isotropic, but the nonlinear 
elasticity is still anisotropic (from the fourth-order terms in the expansion of $ in (2.7) in 
powers of 62 and 63). 

The elastic stress tensor cr= {crij} is expressed as 

a^x = Kei + aijj + fi2 sin(27re2)/27r, 
ayy = Kei + aijj - ^2 sin(27re2)/27r, 

(^xy = cTj,^ = /i3sin(27re3)/27r. (2.12) 

In the linear elasticity, sin(27re2)/27r and sin(27re3)/27r are replaced by 62 and 63, respectively. 
Notice the relation, 

V- a= -^F, (2.13) 
where tp is fixed in the functional derivative 6F/6u. 



The mechanical equihbrium condition V- cr= is equivalent to the extremum condition 
6F/6u = 0. In the coherent case this condition may be assumed even in dynamics. In fact, 
using this condition in the linear elasticity, the elastic field has been expressed in terms of 
ip in the previous theories (see the appendix)i*2jMi^. We then find the following, (i) The 
typical strain around domains is given bj-*^ 

eo = «Ac/Lo, (2.14) 

where Ac = Aip/2 is the composition difference between the two phases and 

Lo = K + fi2o (2.15) 

is the longitudinal elastic modulus. This strain needs to be small (cq < 1/4 approximately) 
as long as the system stays in the coherent regime, (ii) As will be shown in the the appendix, 
in the limit of weak cubic elasticity and weak elastic inhomogeneity, one-phase states become 
linearly unstable for k^lT — Tq + T{ip)'^/2]/vo < / Lq. At the critical composition {ifj) = 
this condition becomes T < Tg with 

% = To + voa'/LokB. (2.16) 

(iii) Furthermore, (A. 5) suggests that the typical domain size in steady pinned states is a 
decreasing function of the quench depth Tg — T. 



III. DYNAMIC EQUATIONS 

In the incoherent case the mechanical equilibrium does not hold around dislocation cores 
when dislocations are created and when they are movingi^. We thus need to set up the 
dynamic equation for the elastic displacement u. In this paper the lattice velocity v = du/dt 
obeys the momentum equation^i, 

p-Q^ = r]oV^v + V- a . (3.1) 

The mass density p and the shear viscosity ?7o are constants. We neglect the bulk viscosity 
term (oc V ■ v) in (3.1) for simplicity^^. In our model the sound waves relax owing to this 
viscous damping and the mechanical equilibrium V- cr= is rapidly attained unless r^o is 
very small. Note that the nonlinear terms in (3.1) are only those in cTjj in (2.12). 



The composition obeys the diffusive equation, 

I^V.AWV- (3.) 
The Icinetic coefficient depends on ip a^^^ 

A(7/;) = Ao(1-V^')=4AoC^cb, (3.3) 
where Aq is a constant. Here u is fixed in the chemical potential difference 5F/5ip, so 

CV V + «ei 



S£ _ kB 

5ip vq 



+ 0[1 - cos(27re2)] + 0[1 - cos(27re3)]. (3.4) 

The last two terms arise from the elastic inhomogeneity. If \{'4') is of the form of (3.3), the 
diffusion equation dcAldt = DqV'^ca {803/ dt = DqV^cb) follows in the dilute limit — >• 
{cb 0) with 

Dq = XokBTvQ^, (3.5) 

where the coupling to the elastic field becomes negligible. In usual solid mixtures the 
diffusion is very slow and vacancies are in many cases crucial for a microscopic description 
of diffusion^Si. Effects of such point defects are not treated in the present theory. 

The total free energy Ftot = F + J drpv'^/2 including the kinetic energy then changes in 
time as 



d 
~dt 



Ftot = - f dr\Y,r]o{ViVjY + A(^) 



5^; 



(3.6) 



Here the surface integrals have been omitted, which vanish if the boundaries are fixed and 
there is no flux of the atoms from outside (and also if the periodic boundary condition is 
imposed in simulations). The above time-derivative is non-positive-definite. As a result, the 
equilibrium is attained when v = V- (t= and SF/dip =const. 

If the lattice is deformed significantly, we should add the convective term —Vlipv) on 
the right hand side of (3.2), treating (3.2) as the equation in the Euler description. If its 
presence is assumed, another term of the form, —ilj'VSF/Sip, becomes also needed on the 
right hand side of (3.1). With these two terms we again have dFtot/dt ^0. These two terms 
are well-known in critical dynamics of fluids^. However, in our solid case, the magnitude 
of the displacement Au = Jq dt'v{r,t') remains small and these two terms give rise to no 
essential differences in our results at not large applied strains^l. 



Also note that the dynamic equations (3.2) and (3.3) may be treated as Langevin equa- 
tions with addition of the random noise terms related to the kinetic coefficients A (■?/') and 
?7o via the fluctuation-dissipation relations^. In this paper, however, we neglect the random 
noise, because the thermal energy k-^T will be assumed to be much smaller than the typical 
energy of elastic deformations. 

IV. NUMERICAL RESULTS 
A. Method 

We integrated (3.1) and (3.2) in two dimensions on a 256 x 256 square lattice. The 
mesh size Ax was set equal to the lattice constant a in the reference state with it = 0, so 
the system length is Lq = 256a. The vectors, u and v, are deflned at the lattice points 
{n,m), while the strains, the tensors, and the composition are deflned on the middle points 
(n+1/2, m+1/2). These are needed to realize well-deflned microscopic slips in our numerical 
schemei^. The periodic boundary condition was imposed except the simulation of applying 
uniaxial deformation (Figs. 13 and 14). Because the time scale of u is shorter than that of 
ip, we integrated (3.1) using an implicit Crank-Nicolson method. Space and time will be 
measured in units of a and 



respectively, where /i2o is deflned by (2.10) and {1^20/ pY^'^ is the transverse sound velocity 
propagating in the [11] direction. The free energies and the free energy densities are measured 
in units of ^20(1^ and ^20, respectively. For simplicity, the scaled time r^^t, position vector 
a~^r, and displacement vector a~^u will be written as t, r, and u, respectively, in the same 
notation. 

In this paper we set K/jj,2Q = 4.5, Q;//i2o = 0.6, and k-QTQ/vQjj,2Q = 0.05, where Tq is the 
mean-fleld critical temperature in (2.4). Since T Tq hereafter, the elastic energy to create 
a single slip (~ ^20^0) is much larger than /cb^ in our simulations. Furthermore, we assume 
weak cubic elastic anisotropy with /i3o//^2o = 1-1 and moderate elastic inhomogeneity with 

1^21 = /^31 = 0.6/i20- 

The dimensionless kinetic coefficients are given by 



To = {p/iJ2oY^'^a, 



(4.1) 



Aq — AoTo/i20'^ 
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^0 = '7o/t-o/^2o- 



(4.2) 



We set X*o = 0.001 and r]* = 0.1. Then, 

A*//i* = Dop/vo ~ 10-^ (4.3) 

Since the relaxation rate of a sound with wave number k is r]ok'^/p, the time scale of i/j 
becomes longer than that of the elastic field by two orders of magnitude. In real solid alloys, 
these two time scales are much more distinctly separated, probably except for hydrogen- 
metal systems where the protons diffuse quickly^. 

In homogeneous one-phase states we have 62 = 63 = and ei = —aip/K. Here well-known 
is a parameter 77 = Ida/dipl/a representing the strength of the composition-dependence of 
the lattice constant a in a mixture^. In our case we have 1] = a/2K = 0.067 and the spinodal 
temperature Tg in (2.16) becomes 2.31. 

B. Slips and composition changes 

Edge dislocations appear in the form of slips or dipole pairsi^, because a single isolated 
dislocation requires a very large elastic energy. Slips are thus fundamental units of plastic 
deformations. In Fig. 2 we show the displacement and the composition around typical slips 
in a one-phase steady state with length lOv^a in the upper plate and 10a in the lower 
plate. Here we initially prepared a slip given by the linear elasticity theory^ at the critical 
composition {{tp) = 0) and let u and ip relax until the steady state was achieved. The 
temperature was kept at T/Tq = 2.5 and no phase separation occurred. As in the previous 
simulation a^^i^^ , we can see Cottrell atmospheres around the dislocation cores. The maxi- 
mum and minimum of ip at the lattice points close to the dislocation cores are of order ±0.6. 
Cottrell's result is obtained as follows: Let T be much higher than Tq and / Lq and the 
gradient term be neglected; then, the condition 6F/6ip =const. yields c^/(l — ca) =const. 
exp{—U/k^T), where U = foaci— . In our case the maximum of lUl/k^T at the lattice 
points is of order 1 and the accumulation is not very strong. 

As a next step, starting with the configuration in Fig.2, we lowered the temperature to 
T/Tq = 2 to induce spinodal decomposition. Subsequently the Cottrell atmospheres grew 
into domains and the dislocation cores stayed at the interface regions. The domain size 
attained finally was of order 50a. Fig. 3 illustrates the displacement and the composition in 
the final steady state, where the maximum and minimum of ip are about ±0.9. Leonard and 



Desaii^ obtained similar composition profiles in spinodal decomposition, where the elastic 
field of dislocations (given by the linear elasticity theory) was fixed in space and time. 

Mathematically, slips in steady states satisfy 6F/6u = and 6F/6ip =const. Without 
externally applied strains, they are metastable owing to the Peierls potential energy arising 
from the discreteness of the lattice structure^^. Although not discussed in this paper, slips 
become unstable against expansion or shrinkage with increasing applied strain. 

C. Dislocation formation around a hard domain 

Fig. 4 shows a single large hard (A-rich) domain at the center in the coherent condition at 
shallow quenching T/Tq = 2 after a long equilibration time. Here is about 0.7 inside the 
domain and about —0.7 outside it. Its shape slightly deviates from sphericity owing to the 
weak cubic anisotropy assumed in this paper. We next performed a second deeper quenching 
to T/To = 1. Subsequent diffusional adjustment of the composition proceeded very slowly, 
but a discontinuity of the order parameter Atp about 1.8 was established relatively rapidly 
across the interface^. As a result, at a time about 1000 after the second quenching, the 
maximum of |e2| reached 1/4, the value at the stability limit, in the interface region (see 
the sentences below (2.7)). We then observed formation of dislocations and generation of 
sound waves emitted from the dislocations. The upper panel of Fig. 5 shows the coherent 
elastic displacement Wcoh just before the dislocation formation, while the lower panel shows 
the subsequent additional incoherent change 6u = u — ttcoh after a time interval of 1000. 
The free energy F in the state in the lower panel is smaller than that in the upper panel 
by 152.9 in units of ^200^ ■ More details are as folUows: (i) two pairs of dislocation dipoles 
(four dislocations) appeared simultaneously in a narrow region, (ii) two of them glided 
preferentially into the softer region forming two slips perpendicular to each other, and (iii) 
slips collided in many cases and stopped far from the droplet, resulting in a nearly steady 
elastic deformation. Thus a half of the dislocation cores stayed at the interface and the 
others were distributed around the domain. These three processes took only a short time of 
order 100. 

After the above dislocation formation at a relatively early stage, the composition changed 
very slowly. We show three figures at t = 23000. In Fig. 6 we displays the following strain, 

e={el + elf'\ (4.4) 



which is invariant with respect to the rotational transformation (2.11). The shps make an 
angle of ±7r/4 with respect to the x axis in the regions with large |e2| (in the uniaxially 
deformed regions), while they are parallel to the x oi y axis in the corner regions with large 
Icsf^. We also notice that the dislocation formation took place with the symmetry axis in 
the [11] direction for our special geometry. Fig. 7 gives the free energy density / in (2.3), 
where the peaks outside the domain represent the dislocation cores. In the interface region it 
exhibits a cliff-like structure arising from the gradient term and higher peaks arising from the 
dislocation cores. Fig. 8 shows the order parameter ip, where we can see Cottrell atmospheres 
around the dislocation cores surrounding the domain. The system is still transient and there 
is still a small composition flux through the interface. 

D. Dislocation formation in a soft network 

Next we examine dislocation formation when hard rectangular domains are densely dis- 
tributed and wrapped by a percolated soft network. As in Fig. 9, we prepared such a steady 
domain structure at T/Tq = 2 in the coherent condition. As in the previous simulation a^^i^^ , 
the hard domains (in gray) are elastically isotropic, while the soft network (in white) is 
mostly uniaxially stretched. That is, in the soft stripes between the two adjacent hard do- 
mains, we obtain 62 ~ 0.2 in the horizontal stripes and 62 ~ —0.2 in the vertical stripes. We 
then quenched T to T/Tq = 1 to induce the composition readjustment. Fig. 10 displays the 
resultant time evolution of the total free energy F = J drf and the snapshots of e in (4.4) at 
the points A, ■ ■ ■ , and E. It demonstrates that F mainly decreases due to the composition 
change but sometimes due to appearance and gliding of slips in the soft stripes. Note that 
the overall composition adjustment occurs slowly on the time scale of R^/Dq = 10^ — 10^ 
where R is the domain size. In Fig. 11 we show the displacement u within the square window 
in B, C and D, respectively, while in Fig. 12 the bird views of the free energy density / the 
square window are given at t = and 4475 after the second quench. Fig. 12 clearly illustrates 
appearance of the peaks representing the dislocation cores. 



E. Uniaxial stretching in two-phase states 



Finally we apply a constant uniaxial deformation to initially coherent states with {ip) = 
to induce plastic flow. That is, we set = Uy = at the bottom {y = 0) and = —Uy = 
eLo/2 at the top {y = Lq). The applied strain rate was fixed at e = 10""^, so e = et with t 
being the time after application of the deformation. In Fig. 13 we plot the average normal 
stress A''i vs the applied strain e for T/Tq = 3, 2.4, and 2 (upper plate), where 

^^1 = {(^xx - (yyy) = -(/i2sin(27re2)), (4.5) 

TT 

where (■ ■ ■) denotes taking the spatial average. The snapshots of e in (4.4) are also given at 
the points a, b, and c (lower plates). For T/Tq = 3 the system is in a homogeneous one-phase 
state and random numbers with variance 0.01 were assigned to ip at the lattice points at 
t = 0. In the initial state at T/Tq = 2.4 the maximum and minimum of ip and 62 are ±0.32 
and ±0.05, respectively. At T/Tq = 2 these numbers are magnified to ±0.75 and ±0.20. 
All the initial states are coherent without dislocations. For T/Tq = 3 the elastic instability 
occurs at e = 1/4 resulting in a fine mesh of slips as in the lower left plate. For T/Tq = 2.4 
the onset point of the slip formation is decreased to e = 0.17. For T/Tq = 3 the onset 
is very early at 0.015, the stress-strain relation exhibits zig-zag behavior upon appearance 
of slips, and the stress continue to increase on the average (up to the upper bound of e 
given by 0.35 in the simulation). Fig. 14 consists of snapshots of e in (4.4) and the shear 
deformation energy density $ in (2.7) in units of fi2o- We can see quadratic appearance of 
dislocations at the center of the uniaxially stretched stripes at e = 0.05 (top plate), gliding 
of the dislocations and pinning at the interfaces at e = 0.1 (top and middle plates), and 
thickening of the slips into shear bands at e = 0.2 (bottom plate). 

We mention a creep experiment in the presence of high volume fractions of 7' precipitates^, 
where softer disordered 7 regions were observed to be filled with dislocation networks after 
large deformations. 

V. SUMMARY AND CONCLUDING REMARKS 

In summary, we have presented a coarse-grained phase field model of plastic deformations 
in two-phase alloys. Though our simulations have been performed in two dimensions, a 
number of insights into the very complex processes of plasticity have been gained. We 



mention them and give some remarks. 

(i) Performing a two-step quench, we have numerically examined dislocation formation 
around the interface regions, which occur spontaneously in deeply quenched phase separa- 
tion. Experimentally^'', dislocation formation has been observed around growing 7'(Al3Sc) 
precipitates at low volume fractions when the radii exceeded a threshold about 20nm. 
Such spontaneous dislocation formation with domain growth has not yet been studied 
theoretically. 

(ii) We have found that dislocations glide preferentially into the softer regions with smaller 
shear moduli and tend to be trapped in the interface regions in agreement with a number 
of observations^. Theoretically, the composition-dependence of the elastic moduli (elastic 
inhomogeneity) is a crucial ingredient to explain the experiments. 

(iii) We have applied uniaxial strain to create multiple slips in two-phase alloys which were 
initially in the coherent condition. The dislocation formation starts in the mostly stretched 
middle points of the soft stripes. A stress-strain curve in Fig. 13 at deep quenching is very 
different from the curves in one-phase states. In real two-phase alloys, a similar monotonic 
increase of the stress without overshoot has been observed, but a considerable amount of 
defects should preexist in such experiments particularly in work-hardened samplea^i^^^. 

This work is a first theoretical step to understand complex phenomena of incoherency in 
solids. Finally, we mention two future problems which could be studied numerically in our 
scheme. 

(i) The composition has been taken as a single order parameter. Extension of our theory 
is needed to more general phase separation processes involving an order-disorder phase 
transition^iiJ^ and to diffusionless (Martensitic) structural phase transitions^ii. 

(ii) Dislocations move under applied strain. The motion is complicated when they are 
coupled with an order parameter and when the time scale of the order parameter is slo w^^i^^ . 
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Appendix 



Here we assume weak elastic anisotropy and weak elastic inhomogeneity in the coherent 
condition in two dimensions, supposing shallow quenching. Then we may eliminate the 
elastic field in terms of ip using the mechanical equilibrium condition V • a = in the linear 
elasticity. We consider the space integral of the last two terms in the free energy density 
in (2.3): AF = / drlaeitp + /ei]. We assume that \fi2i\ and \fi3i\ are much smaller than 
Lq = K + /i2o and that C,a = 2(/U2o/ A^so — 1) is small. Then AF may be rewritten as^i^Siii 

,9 1 



AF = I dr 

+ 



2L 







(A.l) 



J dr g2^\{Vl - Vl)w\^ + g3i^\V,Vyw\ 

where w is obtained from the Laplace equation, 

V^w = t/j- {tjj), (A.2) 

with {^lJ) being the average order parameter. In the first line of (A.l) the bilinear terms are 
written with 

r,ub = -(2aV^o)/^2oea. (A.3) 

The term proportional to r^uh gives rise to anisotropic domains^i,. The second line consists 
of the third-order terms with 

g2 = i^2iay2Ll gs = 2/i3iaV^o- (A.4) 

The third-order terms are known to give rise to pinning of domain growth (and some frus- 
tration effects when g2 and g^ have different signs )^>^. 

In our simulations we set = 2(1/1.1 — 1) = —0.18 and Tcub — 0.0043/i205 so the domains 
tend to become square or rectangular with interfaces parallel to the x or y axis. Furthermore, 
we set g2 = (73/4 = 0.0035/X2o- For f^2i ~ the typical domains in pinned two-phase states 
-Re is given by^ 

i?E~7/[/i2i(Ac)3], (A.5) 
where 7 is the surface tension and Ac is the composition difference between the two phases. 



Thus Re decreases as the quenching becomes deeper. 
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FIG. 1: Normalized shear deformation energy $(e2,e3)//xo for the case /i2 = A*3 = The 
elastically stable regions are meshed with solid lines on the surface. 
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FIG. 2: Displacement vector for a slip (dislocation pair) making an angle of 37r/4 (upper plate) and 
(lower plate) with respect to the x (horizontal) axis in a one-phase steady state at T/Tq = 2.5. 
The arrows are from the initial position in a perfect crystal to the deformed position. The degree 
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FIG. 4: Hard domain in a steady state at T/Tq = 2. The displacement in the square region will 
be displayed in Fig.5. 
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FIG. 5: Upper plate: Coherent elastic displacement ttcoh just before birth of dislocations at deep 
quenching at T/Tq = 1. Right: Incoherent elastic displacement after appearance of dislocations. 




FIG. 6: Snapshot of e in (4.4) after dislocation formation, which is zero within the hard domain 
and nonvanishing outside. The shp hnes end at dislocations. 



FIG. 7: Snapshot of the free energy density / in (2.3). The peaks are located near the dislocation 

cores and the cliff-like structure represents the interface free energy density. 




FIG. 8: Snapshot of the order parameter around a hard domain in the incoherent case obtained 
after a two-step quench. The peak structure at the interface arises because the system away from 
the interface is still in a transient state. The peaks around the dislocation cores in the outer soft 
region represent Cottrell atmospheres (but the minima paired are not seen in the figure). 



FIG. 9: Domain structure obtained at a shall quench T/Tq = 2 in the coherent condition. 
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FIG. 10: Relaxation of the total free energy F in units of 1^20^0 = 20A;b2o after a two-stepnch 
from T/Tq = 2 to 1 with the initial configuration in Fig. 9. It mostly relaxes due to the gradual 
composition adjustment, but it sometimes relaxed due to dislocation formation as enlarged in the 
inset. Snapshots of e at the points A ~ E are given in the lower plates. 
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FIG. 11: Elastic displacement u in the marked regions B,G, and D. 




FIG. 12: Bird views of the free energy density / in units of /X20 = 20A;b?oVq ^ at t = and 4475 
after the two-step quench in Fig. 10. 
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FIG. 13: Stress-strain curves after application of uniaxial stretching e = et with e = 10~^ for 
r/To = 3,2.4, and 2. There is no dislocation at t = 0. Snapshots of e in (4.4) at points a,b, and c 
are given below, which represent slip patterns in plastic flow. 




FIG. 14: Snapshots of e in (4.4) and the shear deformation energy density $ in (2.7) at T/Tq = 2 
for € = 0.05,0.1, and 0.2. 



